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It is shown that the complex d xy component is generated in d-wave superconductor in the magnetic 
field. As one enters superconducting state at finite field the normal to superconducting transition 
occurs into bulk d x 2_ y 2 + id xy state . The driving force for the transition is the linear coupling 
between magnetic field and non zero magnetization of the d x 2_ y 2 + id xy condensate. The exter- 
nal magnetic field violates parity and time reversal symmetries and the nodal quasiparticle states 
respond by generating the id xy component of the order parameter, with the magnitude estimated 
to be on the order of few Kelvin. Parity (P) and time reversal (T) symmetries are violated in this 
state. 



PACS numbers: 74.20.-z, 74.25.Dw 

Symmetry of the order parameter and existence of the 
gap nodes in high-T c superconductors has been one of 
the main experimental questions addressed for the last 
few years. By now the majority of the data support the 
d-wave symmetry of the order parameter |l]-Q. However, 
recent experiments on thermal transport in BiSrCuO su- 
perconductor reveal a number of anomalies, prompting 
the suggestions that the secondary superconducting or- 
der parameter is developed in the external magnetic field, 
thus lowering the symmetry of the initial d x 2_ y 2 to order 
parameter that contains more than one distinct compo- 
nents §,§. 

Based on the data it was suggested that first or- 
der phase transition from original ci ;c 2_j / 2-wave (d) to 
d x 2_ y 2 +id xy (d + id') state occurs abruptly in the mag- 
netic field. Theoretical approach, describing this transi- 
tion as an abrupt bulk generation of the secondary id' 
component upon increased field was proposed It 
was also found that the vortices play an important role 
in the observed anomaly |j 10 1. 

Motivated by these experiments I present here the al- 
ternative approach to the field induced d + id' transition 
in high-r c superconductors. I focus on the high field 
region H ~ H C 2 where vortices are so dense that they 
eventually destroy superconducting state. In this regime 
I still find that the d-wave state can be "distorted" by 
the external field, producing d+id' state in the bulk with 
the intrinsic orbital moment, Figs. (1,2). 

D-wave superconductor has low energy quasiparticles 
in the nodes of the d-wave gap. These low lying states 
have a vanishingly small gap and hence can easily re- 
spond to the external perturbations |ll]] . This fact opens 
up the possibility to generate the second component of 
the order parameter, orthogonal to the initial d-wave 
state, from these "normal" quasiparticles at the nodes. 
This softness of the d-wave state to the secondary com- 
ponent generation in the presence of perturbations, such 
as scattering on the surface |l^,|l3| or the magnetic im- 
purity scattering [ p^ , |l5| is the reflection of the marginal 
stability of d-wave superconductors flq] . 
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FIG. 1. Fermi surface (FS), taken to by cylindrical, is 
shown with the angular dependent Ao(G), Ai(0) gaps. The 
d + id' state is fully gapped because for the complex order 
parameter the particle gap will depend on sum of squares of 
real (d) and imaginary (d') amplitudes. 

Here I consider the case of external magnetic field and 
argue that the normal to superconducting transition in 
the field occurs into bulk d + id' state: d-wave state is 
marginally stable in the field. Magnetic field H\\z vio- 
lates P and T symmetries. The superconducting state 
responds to the field by generating the id xy component 
so that the total order parameter symmetry in the bulk 
becomes d + id' . 

The physical origin of the instability is the bulk mag- 
netic moment (M z ) in d+id' state. First I present a qual- 
itative argument for id' induction by the external field. 
The relevant interaction is the (M Z )B coupling to the 
magnetic field B\\z, where B is the induction in external 
field H. In the pure phase one can think of d-wave state 
as an equal admixture of the orbital moment L z = ±2 
pairs: 
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Ao(8) = A o cos20 = ^(exp(2i9) + exp(-2z6)) (1) 

Here is the 2D planar angle of the momentum on the 
Fermi surface, A is the magnitude of the d x 2_ y 2 compo- 
nent. Motivated by the layered structure of the cuprates, 
I consider 2D d x 2__ y 2 superconductor. In the presence of 
the external field H the coefficients of the L z — ±2 com- 
ponents will shift linearly in B with opposite signs: 

Ao(e)^^((l + 77i?)exp(2ze) 
+ (1 - rjB) exp(-2*e)) = A (6) + i BA^Q) (2) 

where Ai(0) oc rj sin 20 - is the d xy component and r\ 
is the coupling constant. The relative phase 7r/2 of these 
two order parameters comes out naturally because d + 
id' state has a partially noncompensated orbital moment 
L z — +2 . The pure d-wave state has nodes of the gap 
and generated id xy gap will "seal" the nodes making the 
state fully gapped, as is shown in Fig.l. 

To show how complex id xy component appears, con- 
sider a macroscopic Ginzburg-Landau (GL) functional 
which allows the linear coupling between the original 
d x 2_ y 2 order parameter Ao(0) = Ao cos 20 and the field 
induced d xy component: Ai(0) = A!sin20. The GL 
functional contains the linear coupling term: 

F mt = i^A A{B + h.c. (3) 

where r\ is the macroscopic coupling constant, Eq.(^|). 
This coupling is possible only for d + id' and not for 
d + is symmetry of the order parameter. The question 
I address here is to find the conditions when the second 
component is generated in the presence of the external 
field. I find that: a) instability into Ao + iAi state occurs 
at T C (H) greater than the transition temperature for the 
initial transition T®(H) in the absence of the coupling 
Pint- Upon entering the superconducting state the sym- 
metry of the condensate in the bulk is d + id' from the 
outset, as is indicated in Fig. 2. This is the main result 
of this paper, b) id xy induced component is linear in the 
applied field Ai/Aq oc H. 
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FIG. 2. The (H,T) phase diagram of d-wave superconduc- 
tor is shown. The dashed line corresponds to the normal — > 
superconductor transition in the absence of the A0A1 cou- 
pling: H° 2 (T). The solid line is the new phase boundary 
H C 2(T) in the presence of the AoA* coupling in the free en- 
ergy. The d + id' state has higher transition temperature and 
hence the first state to form below H C 2 line is the d + id' state 
in the bulk. The effect of the induced d + id' state can be seen 
in the field cooled experiment, as indicated. The slope of the 
H c2 (T) is linear at small T since the shift ST C (H), Eq.(Q) is 
quadratic in H. I ignore here the effects of possible vortex 
melted phases. The new H C 2(T) is drawn out of scale and is 
closer to the H^ 2 (T) line. In this theory the H = transition 
occurs into a pure d-wave state. 

I start with the GL free energy of the d-wave supercon- 
ductor in the vicinity of H C 2(T) field. The input fields 
are: Ao - d x 2_ y 2-yva,ye order parameter, Ai - d xy compo- 
nent and external field along z- axis H. I will consider 
the infinite two dimensional case, where all the gradients 
are in the (xy) plane only. Here for simplicity I also con- 
sider the case of d x 2_ y 2 + id xy admixture only, ignoring 
any other components of the order parameter, that can 
be mixed in |Fij] . The reason is that d x 2_ y 2 + id xy is 
the only singlet state on the square lattice where both P 
and T are broken and the condensate has nonzero orbital 
momentum (L z ) ~ iA Al + h.c. Q. 

The GL functional density takes the form: 

F = J dr[^{T - T°)|A (r)| 2 + ^|A (r)| 4 + 

^|(AA (r))(^-Ao(r)*)| 2 + f- + F mt + ^A^)! 2 ] (4) 

where the first three terms describe the instability of the 
pure d x 2_ y 2 - wave superconductor, with T° being the 
critical temperature in the absence of the field. Tensor 
Kij is the gradient tensor, which I do not need specify 
right now M). I used Dj = Vj —i^Aj,j = x,y in gra- 
dient energy, Ai- is the gauge potential. The last term 
describes the positive energy shift for the d x 2_ y 2 state in 
the presence of the induced secondary component d xy . 
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The gradient terms for the induced d xy component are 
ignored since I consider the case where Ai is small and 
the gradient terms will lead to the small corrections [^0| . 
Fi n t in Eq.(||) is discussed in more details below. For 
general discussion of the GL terms see [ plp^ ]. 

From GL functional Eq.{|4j) it follows that the coupling 
to the external field, as far as the secondary id! compo- 
nent is considered, is linear in Fi nt . Whereas the stiffness 
term (last term in Eq.(^)) is quadratic. Therefore the 
linear term dominates at small Ai(r) and yields nonzero 
equilibrium value for Ai(r). 

One can easily see that the only second order gradient 
term that couples Aq and Ai is : 



9a* G = leads to: 



D x A D y Al - D y A D x Al + h.c. 



(5) 



leading to, among other terms, — i^i? AqA* + h.c, where 
commutator [D x ,D y ] = i^B. Since there is no other 
terms that couple Ao and Ai in the gradient terms, I can 
write the coupling term without derivatives as is done in 
Eq.(|. 

To show that the terms in Eqs.(||, |J) are indeed full 
invariants, recall that the square lattice group -D4 has 
irreps: A\ (s-wave), A2 ~ z , B\ g ~ x 2 —y 2 - d x 2_ y 2-w&ve 
representation, B 2g ~ xy - d xy representation, all of them 
one dimensional and one two-dimensional representation: 
E ~ (x,y) - p wave p3| . Consider for example B lg and 
B 2g coupling in Eq.(|5|): the the product of A A| ~ A 2g . 
The direct product of two derivatives, transforming as E 
each, is a sum of all four one dimensional representations: 
E x E = A\ g + A 2g + Big + B 2g . The antisymmetric 
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combination D X D V 



D y D x = 



H ~ A 2g and I find 



that true scalar can be formed by taking the product 
of two A 2g representations AoA* and [D. x ,D y ]. Hence 
the interaction term Eq.(|^) is a scalar and allowed in 
free energy. Eq(^) follows immediately. Alternatively, 
as indicated above, one can see that the Ao + iA\ state 
has a nonzero angular momentum state with expectation 
value of the magnetic moment of the condensate (M z ) ~ 
«AqA*+/i.c.. Therefore there is a linear coupling between 
external field and magnetic moment Fi nt oc —(M Z )B, 
Eq.f). 

Next, I use the fact that for H ~ H c2 the local field 
h(r) is homogeneous and d-wave order parameter takes 
the form Ao(r) = Ao/(r), where f(r) is a function of 
the position of the nodes of the order parameter, e.g. 

f( r ) = Ui( z - z i.) ex P(-\ z i\ 2 / ^ 2 h) where z, = x i +iy l are 
the nodes of the order parameter for the case of pure d- 
wave with isotropic gradient term Kij — KSij and In = 
(c/2eH)i p3j. I will not specify the form of f(r) which 
should come out as a solution for particular choice of the 
lattice of the nodes of the d-wave order parameter in our 
case. Instead I note that the Ai field enters quadratically 
into GL functional and one can integrate over Ai field, 
to obtain the effective theory for Ao. 

Treating Ai and A* as independent variables the min- 
imization of the Gibbs energy G = F — 8a x G = 



Ai(r) - — HA {r) 



•O(|A (r)| 5 



(6) 



hereafter I ignore small difference B — H = —4irM z oc 
0(|Ao(r)| 2 ) near H c2 , where Ao is small. After substi- 
tuting this result in Eq. (|j) the free energy for Ao acquires 
negative contribution SF = — f dr[r] 2 /2aiH 2 \A (r)\ 2 ]. I 
obtain the GL functional Eq.(|) with renormalized Tj?, 
which is now field dependent: 



T°(H) = T r ° + ST C (H), ST C {H) = 



(7) 



a ai 

The standard next step is to solve the GL equations for 
the Aq (r), using the ansatz Ao(r) = Aq f(r) as discussed 
above |2l|,^,Q . One would find the free energy for the 
amplitude of the order parameter Ao (now Ao,Ai are 
homogeneous fields corresponding to the amplitudes of 
d + id 1 order) with the field dependent T C (H): 



F = a(T, H)/2\A \ 2 + 0(T, H)/4\A 



,Ai = — HA 

%(X\ 



(8) 



with a(T,H) = a (T - T C (H) - ST C {H)). The a(H,T) 
and P(H, T) parameters will be determined by the vor- 
tex lattice structure. Here the vortex lattice in not spec- 
ified explicitly and I assume that for given configuration 
one finds the solution. I write the general expression 
for thus determined a(H, T) oc T - T C (H) in the ab- 
sence of the AoA^ coupling. It follows from Eq.(||) that 
this coupling will shift amplitude instability to higher 
fields/temperatures T C {H) — > T C (H)+ST C (H) as is shown 
in Fig.©. 

To estimate the amplitude of the induced Ai one can 
use the BCS theory for the coefficients in the GL en- 
ergy Eq. (f|]|) . These estimates are useful as the order of 
magnitude estimates at best. The coefficient r\ in Fi nt 
(in units of energy) describes the (M)H interaction with 
the orbital moment of the condensate and I estimate 
r] = hbNq , where \ib — -j^ c is the Bohr magneton, 
Ao = 1/eV per unit cell- is the density of states at the 
Fermi surface. a± was shown to be a± ~ No/2 [ fL5| . 
Using H c2 ~ 0o/27r£ 2 , with £ = 20 A as the low tem- 
perature coherence length of superconductor, one gets 
77 ~ No(H/ H c2 ){a/£) 2 , where a = 3.8A is the atomic 
length in the Cu-0 plane. From these estimates and 
Eq.(g) it follows that: 

This makes the amplitude of the induced component 
|Ai| to be on the order of few K° . From Eq.(^J) the 
shift of T c turns out to be much smaller: 5T C /T® ~ 
{H/H c2 )(a/tf ~1CT 3 . 

Eqs.(^|j7[) prove the points that a) the transition into 
d + id' state indeed occurs at higher T c if the Ao A^ cou- 
pling is present and b) the field induced A1/A0 xy com- 
ponent is linearly proportional to the field. Whether the 



H/H c2 )(a/0 2 - 10" 



(9) 
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existence of the induced id! component can be observed 
in experiments, such as penetration depth X(H,T), is 
an interesting question which is currently under investi- 
gation. It is possible that the secondary component will 
cause changes in the vortex lattice symmetry at low tem- 
peratures. 

It is useful to note the relation of these results to 
the previous work. Laughlin argued for the coupling 
term similar to F int , Eq(||), although the focus of his 
work was on the different part of the phase diagram 
H c i <C H <C H c2 0. It was predicted that the tran- 
sition into d + id' state to be of the first kind with the 
induced component Ai oc H?. In another approach, 
Ramakrishnan recently suggested the alternative mech- 
anism to generate the d + id' state in the bulk ||. His 
approach is based on the finite id' component generated 
near the vortex cores. As one increases the vortex density 
the overlapped id' patches eventually produce the bulk 
id' order. The starting point there was the single vortex 
solution which is inapplicable at H ~ H C 2- The rela- 
tion between phases discussed in |?],|| and present work 
is unclear at a moment. One possibility is that there is 
a crossover or phase transition that separates high field 
d + id' state, discussed here, from possible d + id' state 
at lower H <C H C 2 fields. 

In conclusion, I argue that the normal-to- 
superconducting transition of the d-wave superconductor 
occurs into d + id' state with the field induced secondary 
component id' with Ai/Ao oc H. Even if the interactions 
are repulsive in the xy channel this field induced compo- 
nent is present at finite fields H ~ H C 2- The driving 
force for the id' induction is the linear coupling between 
orbital moment of the d + id' condensate with the ex- 
ternal magnetic field. Penetration depth A and vortex 
lattice structure could be sensitive to the presence of id' 
gap and might be used to detect it. 

I am grateful to D.H. Lee, V.P. Mineev and M. Zhit- 
omirsky for the useful discussions. This work was sup- 
ported by the US DOE. 
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